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We analyse the spin wave approximation for the 2D-XY model, directly in reciprocal space. In 
this limit the model is diagonal and the normal modes are statistically independent. Despite this 
simplicity non-trivial critical properties are observed and exploited. We confirm that the observed 
asymmetry for the probability density function for order parameter fluctuations comes from the 
divergence of the mode amplitudes across the Brillouin zone. We show that the asymmetry is 
a many body effect despite the importance played by the zone centre. The precise form of the 
function is dependent on the details of the Gibbs measure giving weight to the idea that an effective 
Gibbs measure should exist in non-equilibrium systems if a similar distribution is observed. 



I. INTRODUCTION 

Order parameter fluctuations in critical systems are known to be non-Gaussian ]lj and the form and universal nature 
of the fluctuations has been studied from a field theoretic and renormalization group point of view j^-|5| . More recently 
a surprising and more pragmatic motivation has come from the experimental observation that the probability density 
for fluctuations of global measures in correlated, but non-equilibrium systems, such as energy injected into a closed 
turbulent flow, or global measures of self-organized criticality can be of the same form as those for order parameter 
fluctuations in a critical system . In this context we have made a series of studies of order parameter fluctuations 
in the low temperature phase of the 2D-XY model |l(] jl2| . This is arguably the simplest system in which one can 
study critical phenomena as, in zero field a continuous line of critical points stretches from the Kosterlitz-Thouless- 
Berezinskii phase transition [ jl3| , Tkt, to zero temperature. Critical behaviour therefore occurs even at temperatures 
where a harmonic approximation to the starting Hamiltonian is valid. Further, it follows from a renormalization 
group analysis |l4j that the behaviour at all temperatures below Tkt is perfectly captured at large length scales by 
a harmonic Hamiltonian with effective coupling. As a result, we have been able to calculate the probability density 
function P(m) for all temperatures below which vortices make a measurable contribution, without either invoking the 
scaling hypothesis or using renormalization group theory jll],[l2j . 

In this paper we take further advantage of the simplicity of the model, making a study of P(m) directly in Fourier 
space. The harmonic Hamiltonian is diagonalized in reciprocal space and the problem becomes one of a set of 
statistically independent normal modes. One can thus study the contribution to P(m) from different parts of the 
Brillouin zone. In section II we show that it is essentially the soft modes near the zone centre that are responsible for 
the asymmetric non-Gaussian form for P(m). However, we show further that a complete and detailed description of 
the P(m) requires many modes and the problem remains intrinsically a many body one. 

In section III we discuss our results in relation to extremal statistics. We have previously fitted P(m) to a function 
reminiscent of Gumbel's first asymptote from extremal statistics. With this description in terms of normal modes 
we are in a position to study the contribution from the mode with the largest amplitude in each configuration. We 
confirm that it is not sufficient to keep the largest one of the statistically independent modes and conclude that, if 
extremal statistics are relevant for the description of the global measure, then they must be applied to composite 
correlated structures and not to the most fundamental excitations of the problem. 

In section IV we study the effect of changing the form of the microscopic generator for the normal modes. We show 
that changing from the Gibbs measure to an arbitrary form does change the final outcome for the global fluctuations. 
We generalize the Gibbs measure analytically using a x 2 (V) distribution for the contribution from each mode. This 
function gives the probability for a composite variable constituting a sum of v statistically independent degrees of 
freedom. We show that the global measure becomes Gaussian as v becomes large and give a physical interpretation 
for this. Some conclusions are presented in section V. 
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II. DESCRIPTION OF THE TWO-DIMENSIONAL XY MODEL IN FOURIER SPACE 



A. Framework of the study 



In the 2D — XY, or plane rotator model classical spins of unit length, confined to a plane in spin-configurational 
space interact via the Hamiltonian 

H = -J & ■ $3 =- J cos ^ ~ e i)i w 

<i,j> <i,3> 

where the orientation of the spin vector is give by the angle 6i- J is the exchange constant, the sum is over nearest- 
neighbours and here we consider spins on sites of a square lattice with periodic boundaries and lattice constant of 
unit length. 

At low temperature the harmonic, or spin wave Hamiltonian can be written 
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where <f> q — Re l/s/N'^2 i 9iexp(— iq.fi) and j q = 4 — 2cos(q I ) — 2cos(gj,). Here (j> q takes on values between ±oo 

as we have neglected the periodicity of the original spin variables. The sum is over the N — 1 modes of the Brillouin 
zone (BZ) neglecting the Goldstone mode at q = 0. The latter is a purely diffusive mode contributing zero in energy. 
As discussed below we find that we can replace j q = q 2 without introducing an error in the calculated distribution 
function. 

In principle this Hamiltonian correctly describes the critical behaviour throughout the low temperature phase as 
bound vortex pairs are renormalized away at very large length scale. In practice this is true for the finite system sizes 
studied here only as long as no vortex pairs are excited. However, as vortices only appear in a relatively small range of 
temperature close to the KTB transition JTs) there remains a large range of temperature over which the behaviour is 
correctly described by a Hamiltonian of this form. Within this regime non-linear but analytic terms in the expansion 
of cos(6>i — 9j) are accounted for perfectly by an effective coupling constant QJllJ. 

Thus, at low temperature, the 2D XY model can be viewed as a polydisperse perfect gas of particles of "mass" 
Jj q . The probability function for each mode is therefore a Maxwell-Boltzmann distribution 



P{<Pq) = y eX P( J}""*) ® 

depending explicitly on the vector q. 

The aim of this article is to show that this dispersion of mass, associated with a Gibbs measure is at the heart of 
the mechanism of criticality and that it leads directly to the asymmetry of the probability density function (PDF) 
for the magnetic order parameter. Given the separable nature of the problem we are able to do this by performing 
numerical simulation and analytic work directly in Fourier space, without referring to the spins in real space. The 
order parameter m for a single configuration can be defined 

1 N 

i=l 

where 

is the instantaneous magnetization direction. With this definition, m is a scalar quantity giving the length of the 
magnetization vector independently of its direction in the spin plane. 

In a previous paper (ll), we have calculated the PDF for a large, but finite size system, as a function of the reduced 
variable fi = (m — (m))/a, where (to) and a are the mean and the standard deviation of the distribution: 
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dx 



Pl(p) = I ^e^(x) (6) 
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where Pl(jj) = <xP(m), G(q) = 7" 1 and gk — J2bz G(q) h /N k . The sum over q and the integral over x, in (||) 
are performed numerically and Pl(p-) approaches a (thermodynamic ) limit function rapidly for N > 100 spins. 
Surprisingly this function, shown in Figure 1 and discussed in detail below, is independent of temperature and 
therefore of critical exponent rj = T/2nJ along the line of critical points. 

Expanding (||) in small angles, keeping only the leading terms and Fourier transforming, one can write 

BZ BZ 

where m q = -^fq 2 has the normalized distribution 

P W = ^M m -l/2 e -^, . (8) 

While this expansion is strictly valid to first order in T and for finite systems only, we have previously shown that 
the PDF for the fluctuations in m, given by (Q) is identical to that for m given by (Q), even outside the range of T 
for which the expansion is valid. For this linear development of the order parameter, the critical exponent, rj = 0. 
There are a number of subtil points associated with this result; most notably, the violation of the hyperscaling relation 
between critical exponents if definition (^) is used. These points will not be discussed further here; rather, we refer the 
reader to references JlT| , ^2| ,p[ where they are discussed in detail. However the important point for this paper is that 
the magnetic fluctuations of the entire low temperature phase can be characterized by the harmonic Hamiltonian (|^) 
and the linearized order parameter (J?]) and the problem can be analysed entirely in Fourier space, with statistically 
independent harmonic degrees of freedom. 



B. Simulations in Fourier's space 

We first test these ideas by comparing the PDF, as generated by @, against Monte Carlo simulation performed 
directly in q-space. We take the real part of the Fourier transformed spin variable as the independent degree of 
freedom on each point in the Brillouin zone and define the quarter of the zone with q x and q y positive. Points with 
q x 7^ and q y ^ have a 4-fold degeneracy, points on the axis, with q x = or q y = have a 2-fold degeneracy and 
the origin is excluded, giving the total of N — 1 points. For each configuration, amplitude (p q is generated using (||). 
The results of the simulation are compared in Figure 1 with the theoretical result given by equation (^). Data are 
shown for 10 s configurations. Agreement is quite excellent, showing clearly that the model of statistically independent 
variables and the calculation leading to (||) are consistent and that they correctly describe the fluctuations in the low 
temperature domain. The PDF is asymmetric, with a quasi-exponential tail for fluctuations below the mean, and 
with a much faster double exponential fall off for fluctuations above the mean [O . 
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FIG. 1. The PDF in two dimension (L = 32) at temperature T = 1. The continuous line is the fast Fourier transform, 
equation (ph. The dotted line is for q-space Monte Carlo. 
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C. Relevance of soft modes 



Our aim is now to study the effect on the PDF coming from different parts of the Brillouin Zone and in particular the 
contribution from the soft modes, near the zone centre. We should be able to do this by restricting the summations 
in equation (0) to the desired part of the zone only. We define shells moving out from the zone centre: the first 
contains only the four softest modes. The subsequent shells are squares of length and unit thickness, containing 
(2n + 1) 2 — 1 modes. We have tested the PDF generated from equation (^), with restricted sums over the shells defined 
above, against Monte Carlo data taken for modes in the same shells. We again found essentially perfect agreement 
giving confidence in the theoretical result for both complete and partial sums over the Brillouin zone. In the following 
analysis we need to make no further reference to numerical simulation and can restrict ourselves to the analytic result 
of equation (||) . Some initial steps in this approach were made in |l5| for the PDF of height fluctuations in models of 
interface growth. 

Results are shown in Figure 2 for summations over varying numbers of shells in the Brillouin zone and in Figure 3 
for a calculation including all modes except the four softest modes. One can see that the asymmetry comes explicitly 
from the soft modes near the zone centre. The exponential tail is reproduced to a first approximation with only the 
first four modes, while excluding them from the distribution gives a nearly symmetric distribution. 

The reason for this is the dispersion in the typical amplitudes for modes over the Brillouin zone. For small q values 
7 9 ~ q 2 , which means that the average value, (m q ), is of order 0(T/J) at the zone centre, while (m q ) ~ T/NJ on the 
zone edge. The dispersion in contributions therefore diverges with systems size. The magnetization is the sum over 
variables m q ~ (p 2 , not the ip q themselves and the distribution for the m q , (equation ||), therefor have exponential 
tails. If terms in the sum are not individually negligible, as is the case here, the distribution for the global measure is 
consequently asymmetric. However it is not the case that the the global PDF is reproduced quantitatively by softest 
modes by themselves. More modes are required for a quantitative fit and one must include more than the 4 shells 
shown in Figure 2 for a good fit by eye. 

The contribution to the PDF from any set of N g modes with fixed q is therefore that for an ideal gas of N g identical 
particles; the so called \ 2 distribution 

f(x) ~ a;"'/ 2 - 1 exp (-0x), (9) 
Writing the % 2 distribution in terms of \x one finds 

P L (ji) = M(s - m)^/ 2 " 1 exp (&0u - a)) (10) 

where 



ir (aa\ N 9/z 1 1— < M > , aa a 1 N„T , . 

^=Ur ^nvTT. s = Z . a = 8ir 2 J a = -J (11) 



TJ r(^)' " o V 

The PDF is therefore a convolution over \ 2 distributions with varying particle number and "mass" . Once the 
discrete nature of the lattice is lost the number of particles for each value of q is proportional to the density of 
states, n(q) ~ <7 d_1 . For numerical simplicity we have defined square shells but our findings are compatible with this 
continuum description. 

We study the goodness of fit for n shells, in Figure 4, where we plot the error function 

^"i^-^n (12) 

with N e ff = (2n + l) 2 — 1, for N — 10 4 sites. The further out in the wings one goes, the more shells one has to 
include in order to get a required value of 8. For large /i and N > > N e f j we find 

8( f x,N eff )~l + C(»)-?— (13) 

N eff 

independently on N. The many body nature of the distribution is therefore confirmed by the fact that 8 only goes to 
zero as N e ff diverges. 

We have previously argued [jllj that the critical nature of the 2D-XY model is epitomized by the fact that all length 
scales from the microscopic to the macroscopic and therefore all q values are important. For example, calculation of the 
mean magnetization (to) involves an integral of the form J 2w / L 1/qdq ~ log(L) where both limits of integration come 
into play. One could argue that the above result, with 8 becoming independent of N in the limit N >> N e ff >> 1 
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is not quite in agreement with this discussion as only the zone centre is important in this limit. However, exactly the 
same would be true for a situation where the central limit theorem is satisfied, that is, with zero dispersion in the spin 
wave stiffness of the modes: here one would build a Gaussian distribution by adding the modes from different shells 
and in the above limit, this would be attained with arbitrary accuracy by increasing N e ff independently of the value 
of N. The latter is realised in the presence of a strong magnetic field M] where the dispersion in spin wave stiffness 
is removed by the finite magnetic correlation length. The above analysis is not therefore the most sensitive measure 
of criticality, although it docs clearly show that the global measure is the result of many body contributions, while at 
the same time showing why no error is incurred when replacing j q by q 2 in the various summations. 




FIG. 2. PDF generated from the first shell (4 modes, dot-dashed), 1 shells ( 8 modes, long dashed) 2 shells (24 modes, short 
dashed), full Brillouin zone (solid line). All results are generated from equation (0) for L = 100 
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FIG. 3. PDF generated from all but the first shell (N — 5 modes, short dashed), full Brillouin zone (solid line). All results 
are generated from equation (0) for L = 100 
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FIG. 4. 8(6) = ^g^f i)^ for a given size N and for different N e ff ( different fractions ). From top to bottom on the left of 
the figure, the curves are for 4, 8, 25 and 48 modes. 



D. Comparison with the three dimensional XY model 



Criticality can be seen more clearly by comparing with the magnetic fluctuations in the low temperature phase 
of the three dimensional XY model. Here there is broken symmetry and long range magnetic order, although the 
longitudinal susceptibility remains weakly divergent in the ordered phase and temperature is a dangerously irrelevant 
field near T — [ 17|Jl8[ . As a result the PDF for the magnetic fluctuations is weakly asymmetric jll] as shown in 
Figure 5. Equation (|6|) is easily extended to three dimensions and the same partial summations over the Brillouin 
zone can be performed. Included in Figure 5 are the PDF calculated from the softest modes only and calculated from 
all but the softest modes. As one can see, the form of the curve is well reproduced by the first shell. This should be 
compared with the two-dimensional case where many more modes were required for an equivalent description. 




FIG. 5. PDF for the 3D-XY model for the first shell (6 modes, short dashed), all modes except the first shell (long dashed), 
full Brillouin zone (solid line). All results are generated from equation (^) for L = 100 

This can be seen in more detail by considering the evolution of the asymmetry of the distribution through the 
skewness, 7 Q, 
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as one consecutively removes shells around the zone centre from the summations in g^. That is, we now sum over 
the Brillouin zone except for an ever increasing number of shells around the zone centre. The results are shown in 
Figure 6 for two and three dimensions. Clearly it is necessary to remove fewer shells in three dimensions than in two 
dimensions, before 7 — > and the distribution becomes Gaussian. 




2 4 6 8 10 12 
number of infrared shells that have been removed 



FIG. 6. Evolution of the skewness 7 as a function of removed infrared shells. Study in two (full line) and three (dashed line) 
dimensions 



III. LINK WITH EXTREME STATISTICS 



In |7|,|1 1| , it was shown that the analytic expression for the PDF Pjr fu) is extremely well represented by a function 
reminiscent of Gumbel's first asymptote for the statistics of extremes jl9| 



(15) 



where 



a a b 



1 

s= b 



I 1 d 2 T(a) 
r(o) da 2 

1 

log(o) 



1 dT{a) 



-1 2 



r(o) da 
dT(a)~ 



T{a) da 



(16) 



The function therefore has a single independent parameter a. When a is an integer, equation ( [15] ) gives the PDF for 
the a th smallest value from a set of N random numbers in the limit N — > 00. We have made a comparison between 
the function ip of eqn. (^) and the characteristic function for Pq 



Impcix) = In 



wT(o 



*(o) ln(a) 
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*'"(a) 



120k 5 



b 

* (4) (a) 



2b 2 



*'(a) 



'6b 3 



(17) 



where ^(z) is the digamma function T / (z)/T(z). Expanding ip and equating the first four terms we find an implicit 
function 



W"(g) 

*'(o)3/a 



_ 2 3/2. 



93 



.92 



3/2 ' 



(18) 



Solving for a and hence for the other constants we find a = 1.58, b = 0.93, s — 0.37, u> = 2.16. This does not represent 
the exact solution for Pl : higher order terms in the expansion of the two functions are not the same and their ratio 
diverges slowly from unity with increasing order. 

Although (15) is not the exact solution, mearly a good approximation and the best fit value a ss 7r/2 is non-integer, 
it is still interesting to search for a concrete link with extremal statistics. As we have reduced our model to one of 
statistically independent Gaussian variables the obvious question seems to be: is the PDF of the dominant mode the 
same as that of the global measure? We define a magnetization m = 1 — m max , where m max is the largest of the 
Tdq . We compare, in Figure 7 the PDF for this quantity with that for the full magnetization. The difference between 
the two distributions is not very great, but they are certainly not the same and we confirm what we have seen in the 
previous section, that the correct distribution comes from an ensemble of statistically independent objects and not a 
single mode. The extreme values are well described by the Gumbel function with a — 1 , despite the dispersion in the 
spin wave stiffness. The diverging dispersion means that the largest mode is always one of only a few modes near the 
zone centre and it is by no means clear that Gumbel's asymptote should apply here. 

We can conclude therefore that an extremal description is not applicable to the fundamental excitations in the 
problem. However, one can perhaps apply extremal statistics to the complex and correlated structures that appear in 
real space for critical systems and for correlated systems in general. Indeed, in simulations of the 2D Ising model at 
a temperature T* slightly below the critical temperature, PDFs were found to be very similar to that for the 2D-XY 
model, when calculated for both the magnetization and the magnatization coming uniquely from the largest cluster 
of connected spins @ . There are clearly many open questions here but the above analysis strongly suggests that an 
observation of extremal statistics is a consequence of the presence of composite and correlated objects, rather than 
being a directly observable phenomenon among the most fundamental microscopic excitations that make up these 
objects. 
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FIG. 7. Contribution from the Fourier mode of extreme amplitude (dashed line) and comparison with the PDF generated 
from equation (pl) (solid line) and the Gumbel asymptote, a = 1 (dotted line). 



IV. INFLUENCE OF THE MICROSCOPIC MEASURE ON THE GLOBAL MEASURE 

One might also ask the question: is the global PDF dependent on the existence of a Gibbs measure, or effective Gibbs 
measure for the fundamental excitations that contribute to the global fluctuations? Within the present framework 
this question can be addressed directly by changing the form of the generating function (||), while keeping the same 

dispersion for the standard deviations, a q = «/ (m 2 ,) — (m q ) 2 — yfjj^ - - I n Figure 8 we show the global distribution 

produced from an array of flat generators 

p{m q ) = l/x q , m q < x q , (19) 
p(m q ) = 0, m q > x q , 
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with x q — \/\2a q , which satisfies the above criterion. As can bee seen, the resulting distribution is radically different 
from that for the Gibbs measure. It is prefectly symmetric despite the dispersion in mode amplitudes. The curve is 
almost, but not quite Gaussian. We have not made a detailed analysis of this point, but the deviation from Gaussian 
behaviour suggests that the central limit theorem is still violated even though the distribution is now symmetric. 
From this exercise we conclude that the dispersion in normal modes is not the only criterion for the asymmetric 
distribution we have observed. Our results suggest that further to this, we require microscopic generators that are 
themselves asymmetric. In fact, one can deduce this general result by considering Gaussian generators for the m q : 
here, by definition Pj,(m) must be Gaussain, as the sum of any set of Gaussian variables must itself be a Gaussian 
variable and therefore symmetric. 
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FIG. 8. PDF calculated numerically from the flat generators, equation (|19|- symmetric full line) compared with the PDF from 
equation (^asymmetric full line) and a Gaussian distribution (dashed line). All data is for L — 100. 

One can see analytically how the global PDF varies from the original distribution to a Gaussian distribution by 
considering generators of the form 



Pa (m q ) = ^mJ-'e-^/W. 
Here a > in order to ensure the convergence of p(m q ) and N a 



(20) 

is a normalization constant. We again take 



r(a)m° 

(m q ) = l/j3JNj q so that p a =i/2(mq) accounts for the two dimensional XY model with Gibbs measure. 
We proceed as previously, defining m within a linear approximation only, for simplicity. We define 

9^0 



(21) 



and generate the m q using equation (20). The S function imposing the constraints m = 1 — Yl q m q can be rewritten 

; 1 m— (1— > r 



6(m- (l-^mg)) 



dk ik[ 



(22) 



As modes are independent the conditional probability p(m gi , ...,m qN ) is just the product 

p{m qil ;/;,, 

so that 

P(m) -- 



"qi j •••! "'gjv I 
N 

■ - " "qm ) — Pa { m q ) 
9=1 



dk 

2n' 



TT / dm q p a (m q )t 

q=l J ° 



(23) 



(24) 







After some straightforward algebra, one can now express the PDF Pl (fi) 



J — I 



dx 
27 
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where 



\og(ip a (x)) = (-is 



— Ti aTr log(l 
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which is just a generalization of (Jg) for arbitrary values of a, for which 

(m) = 1 - aT^2^rr, cr = T^/ag^. 



(25) 



(26) 



(27) 



As previously, the expression ( |25| ) can be transformed numerically for large but finite N. This procedure can again 
be tested by comparison with Monte Carlo simulation performed directly on the generators m a , The results, which 
we do not show here, agree as well as they did in the initial problem. In Figure 9 we show Pl{h) for various a values 
calculated from equation ((2^). One can see that increasing a reduces the asymmetry of the global distribution, with it 
becoming Gaussian as a becomes large. This result can be interpreted physically by noting that the generator (20) is 
a x 2 distribution for the variable m q , which is the sum of v = 2a statistically independent positive definite variables. 
One can therefore think of the global PDF as being constructed from v = 2a replicas, or identical systems. As a 
increases, the number of contributions coming from each shell in the Brillouin zone, even those near q = increases. 
The contribution from each shell will become Gaussian and Pl(to) will therefore become Gaussian. 
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FIG. 9. Influence of the microscopic distribution p a {m q ) on Pl(^), calculated result from equation (|2q-full line) generalized 
Gumbel function fitting function (equation ^-dashed line). Data sets for decreasing skewness are for a — 1/2, 2, 5. 

The analysis in terms of the Gumble distribution can also be extended to include this variation in the form of the 
microscopic generators. Repeating the analysis of section [II for the generating function ip a we find 



*'(a) 3 / 2 



53 



1 



.92 



3/2 



(28) 



sb = ln(a) - -0(a) w = ^ 



The knowledge of a again allows one to evaluate the other constants b, s and w 

b = V^(a) -^-^ .'.'^ •«- a " b 

and numerically we find 
• a = 2 : 

a = 5.45802, b = 0.44835, s = 0.21054, w = 96.563 
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• a = 5 : 

a = 13.0845, b = 0.2818, s = 0.1373, w = 194537.9 



Data generated from equation (|T^) is included in Figure 9 for comparison with the exact solution. Agreement is 
very good, indicating once again that the modified Gumbel function is a very useful working tool for data analysis. 
The small deviation between the Gumbel function and equations (^5|), observed along the exponential tail of the 
distribution for a = 1/2 disappears as a increases. 

For large values of a, one can expand ( |28| ) using Stirling expansions for the T function and its derivatives and find 
that 

ago „ 2a 

We actually find that this expression is quite accurate even outside the range of a values for which Stirling's formula 
is strictly valid. 

This analysis is very similar to that recently made for the 2D-XY model in the presence of a magnetic field jl2) . 
In that case the field introduces a correlation length, £, which allows one to think of dividing the system into order 
(L/£) 2 statistically independent parts, just as we have proposed above for a different from 1/2. 



V. CONCLUSION 



In this paper we have exploited the appealing property of the 2D-XY model, that non-trivial behaviour of a corre- 
lated systems appears despite the model being itself diagonalizable at low temperature. The magnetization can hence 
be described in terms of a sum over statistically independent variables. The PDF for order parameter fluctuations can, 
in consequence be written as a convolution over x 2 distributions for statistically independent variables, or particles of 
continuously varying amplitude. We have shown that, as one might expect for a critical system, it is the softest modes 
that influence the distribution most, with their diverging amplitude being directly responsible for the violation of the 
central limit theorem and the observed asymmetry for the PDF of this global quantity. However, we also explicitly 
show that the observed asymmetry is not a single or even a few-body phenomena. Rather, a quantitative reconstruc- 
tion of the PDF requires the sum over a number of modes that diverges in the thermodynamic limit. We have also 
analysed the "extreme value" contribution from the largest mode from each configuration. The resulting distribution 
is consistent with the standard theory of extremal statistics and quantitatively different from the the global PDF. 
Reconstructing the PDF by summing successive shells leading out from the zone centre requires an arbitrary number 
of shells for arbitrary accuracy. 

We show that the detailed form for the distribution depends on the precise form for the microscopic distributions for 
the individual modes. Hence we show that the form of the global PDF depends on a second criterion, the existence of 
a Gibbs measure, and that the nature of the violation of the central limit theorem depends on the form of microscopic 
measure. This discovery is perhaps surprising given our empirical observation that a very similar distribution can be 
observed for complex systems driven far from equilibrium and in particular that for injected power fluctuations in a 
closed turbulent flow at constant Reynolds number || . This extra exigence seems to give weight to the idea that an 
effective temperature should exist in such non-equilibrium, steady state systems. More work in this area is clearly 
required. 
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